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Abstract 

Prime-based ordering which is proved to be admissible, is the encoding of indeterminates in 
power-products with prime numbers and ordering them by using the natural number order. 
Using Eiffel, four versions of Buchberger's improved algorithm for obtaining Grobner Bases 
have been developed: two total degree versions, representing power products as strings and the 
other two as integers based on prime-based ordering. The versions are further distinguished by 
implementing coefficients as 64-bit integers and as multiple-precision integers. By using prime- 
based power product coding, iterative or recursive operations on power products are replaced 
with integer operations. It is found that on a series of example polynomial sets, significant 
reductions in computation time of 30% or more are almost always obtained. 
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1. Introduction 

In 1965, Buchberger developed the method of Grobner bases for solving systems of 
multivariate non-linear polynomials. Since then, computing power has grown, and im- 
proved algorithms developed, but, even now, the method remains impractical for many 
problems. The main approach to improvement has been to develop algorithms which 
avoid unnecessary computation culminating in Faugere's F5 algorithm. On the other 
hand, various methods of representing polynomials have been explored and the impact 
of data structures on algorithm performance evaluated. A central element of the method 
is that it is based on the ordering of power products (known as terms in Faugere's pa- 
persfl]- Orderings that have been used include lexicographic, total degree and variations 
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on these. Prime-based ordering does not appear to have been exploited, and it is the 
purpose of this paper to explore this case. 

The plan of this article is as follows. First, we introduce prime-based ordering as 
the natural ordering of power products imposed by encoding the indeterminates as dis- 
tinct prime numbers. Then, we show that this ordering is admissible. We report our 
implementation of Buchberger's improved algorithm using both total degree ordering 
and prime-based ordering. Experimental measurements show that significant gains are 
achieved by using prime-based ordering. 



2. Prime-based ordering: an admissible ordering based on prime numbers 

Total orderings used in Grobner Basis Algorithms are required to be admissible. That 
is they must satisfy the conditions: 

Wt:t^l:l<t (1) 

s<t^s-u<t-u (2) 

Common admissible total orderings include lexicographical, total degree lexicographi- 
cal and degre e rever se lexicographical. Variations of these and other orderings are possible 
Buchbergerl (|l985l) ]. 



An ordering that, to the best of our knowledge, has not been used in implementations 
of Buchberger's algorithm for Grobner Bases is one based on prime numbers. Given a 
power product, t = ^" , each indeterminate Xi is mapped to a unique prime 

number: 

xi ^ 2 

X2 ^ 3 (3) 



so that t ^ rit ^ N . For example, x^y'^z ^ 2^3^5 — 360. 
We define 

s < t 4^ Us < nt (4) 

In other words, the ordering of s and t is determined by the natural ordering of the 
integers Us and nt- We call the integer nt the prime image of t. 

Now, if all the exponents of s are zero, s = 1 and Us = 1. Similarly, if < ^ 1, at least 
one exponent of t is positive, so that nt > 1. This establishes condition (1). 

Similarly, from (4), s < t ^ < nt ^ ngnu < ntnu =^ s ■ u < t ■ u, establishing con- 
dition (2). 

Hence, the mapping of indeterminates onto a set of prime numbers is an admissible 
ordering. Prime-based order is neither a total degree nor a lexicographical ordering. For 
example, in total degree, x^ > y^] but, in prime-based ordering, a;'^ <-> 8 < 9 <-> y^; so 
x^ < y^. Similarly, x^ < xy in lexicographical ordering; but, .t^ <-> 8 > 6 «-> xy, implying 
x^ > xy in prime-based ordering. 
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3. Implementations of Buchberger's Improved Algorithm 



We have developed four versions of Buchberger's improved algorithm [Buchberger 
(|l985| )]. which generate a reduced Grobner Basis. They are written in the object oriented 
programming language Eiffel so that the algorithm is the same but the implementations 
of power products and coefficients differ in each version. Hence, variations can be readily 
created and verified not only by comparing output but also by using Eiffel's built-in 
preconditions and postconditions. Eiffel allows these conditions to be discarded during 
compilation so that the algorithm can run at full speed. Furthermore, mathematical 
constructs can be represented as structures of interconnected objects. 

The versions differ in the power-product representation and the integer precision used 
for the coefficients. In the first version, monomials are represented as structures consisting 
of a rational coefficient and a power product, t — . The coefficient is repre- 

sented by a pair of 64 bit integers - numerator and denominator - and the power product 
is represented as a string of characters, in expanded form. For example, x^ipz is stored 
as the string "aaabbc" . Power products are compared using total degree in this case. 
This representation was chosen for simplicity, and, because object oriented programming 
methods are used, may be changed easily to any other form of representation, such as 
a vector of powers, [3, 2, 1] or as a list. Operations on power products are implemented 
as iterations or recursions on the underlying data structure. For example, multiplying 
x^y^z by xyz involves merging the strings "aaabbc" and "abc" to yield "aaaabbbcc", 
representing x^y^z'^. 

In the second version, primes are used to represent the indeterminates. Power products 
are ordered by their integer image. This avoids iteration or recursion in the basic opera- 
tions, which reduce to integer operations. For example, multiplication of power products is 
integer multipHcation, as in x^y'^z x xyz ^ 360 x 30 = 10800 = 2''3^4^ <-» x'^y^z'^. Simi- 
larly, the functions of division, lowest common multiple and greatest common divisor also 
reduce to integer operations. The same is true for Boolean operations such as compari- 
son, equality, divisibility and so on. The prime based implementation of power-products 
is the same as the original total degree version, but for replacing the routine bodies with 
the prime based equivalent. 

As very large coefficients are often generated, the Gnu multiple precision library, GMP, 
is used to create the other implementations of the algorithm. This was done by changing 
the coefficient implementation to use the multiple-precision integer type in the GMP 
library. 



4. Validation 

Steps were taken to demonstrate that our version of Buchberger's improved algorithm 
generates Grobner Bases. The necessary and sufficient condition to be met by a Grobner 
Basis, F, is : 

V/i, f2 ^ F :: NormalForm{F, S Polynomial{fi, f2)) — (5) 

where NormalForm{., .) is the normal form or reduction function, and SPol ynomial(., .) 
is leading term cross-reduction function as defined by Definition 6.4 in (iBuchberger 
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(|l985f )]. In the case of a reduced Grobner Basis, in addition to the above, it is also 
necessary that 



V/ e F :: NormalForm{F - {/}, /) = / (6) 

All sets of polynomials generated by the algorithm were tested for these conditions. It 
was also verified that each result set reduced the given polynomial set, confirming that 
the Ideal was unchanged, and that the basis generated by the Grobner Bases package 
supplied with Maple reduced the result, hence providing independent confirmation. 

We found one problem with the implementation of the algorithm as presented in 
Buchberger1(|l985l) ]. This was that the set of pairs of polynomials was not updated after 



a new polynomial was generated by reduction in the Subalgorithm NewBasis{., ., .). To 
the best of our knowledge, this has not been previously pointed out. When this deficiency 
was remedied, results proved to be consistent. 



5. Experimental Results and Discussion 

Our timing tests have been performed using the four implementations as shown in 
Table 1. Examples 1, 2 and 3 are simple manufactured examples, while most of the 



others are take n from [Giovini I (|l991r )] . The examples entitled "Arnold" are taken from 
Arnold I ()2003l) ]. 

Care was taken in the ordering of the indeterminates in the case of "Parametric Curve" 
as it includes the power product x^^. This power product can be encoded as 2^^ when 
using 64-bit integers, but it cannot be encoded as S'^^, which requires 72-bits. 

The main result is that representing the indeterminates by prime numbers and the 
power products by unique integers significantly reduces the computation time in most 
cases. In eight cases, the computations complete successfully when the coefficients are 
encoded as 64-bit integers, and the reduction in time is at least 30%, or a 40% speedup. 
The greatest reduction, ("Gerdt 1"), is 96.8%, or more than a 30x speedup. As there 
are more efficient coding schemes than the one used for the total degree implementation, 
these speedups are on the optimistic side. 

A second result is that the number of polynomials may be different. For example, 
"Gerdt 1" reduces to 56 polynomials when total-degree ordering is used. When prime- 
based ordering is used, only 36 polynomials are generatecPI. However, this is not a general 
rule, as "Gerdt 3" generates 21 polynomials for total-degree and 23 for prime-based 
ordering. 

For some examples, integer overflow problems arise because the coefficients are encoded 
as pairs of 64-bit integers. To counter this, the Gnu Multiple Precision (GMP) library 
is used to support coefficients based on very large integers. Using the library allows 
eleven cases to be compared because they complete successfully with both prime-based 
and string based power products. A twelfth was completed when garbage collection was 
turned off. 



Plex ordering, however, generates 26 polynomials for this case. 
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Table 1. Experimental Results, showing the execution time in ms and 



the number of polynomials in the generated basis 



T'lTTlP f TTm 1 / 
-I- IXXXtJ 1 ±±±o / / 

Number of 
polynomials 


64-bit integer coefficients 


Multiple-precision integer coefficients 


Power Products 


Reduction 


Power Products 


Reduction 
(/oj 


total degree 


prime based 


total degree 


prime based 


Example 1 


13.22/1 


9.19/1 


30.5 


29.92/1 


34.09/1 


12.2 


iLxampie a 


2.93/6 


1.98/6 


32.4 


7.32/6 


6.4/6 


12.6 


Example 3 


9.44/6 


6.35/6 


32.7 


24.82/6 


20.26/6 


18.4 


L/ycuc ^ 


10.73/7 


5.43/7 


49.4 


21.4/7 


15.54/7 


27.4 


V \ C ■) 


1280 5/20 


a 




1) 


142<S!)/24 




VtCKU 1 


.'U-^TOO/^fi 




9().<S 


■■-')(S7()()4/r)() 


1497-")/;-')() 


!)().! 


Liercit A 


56.8/8 


4.91/5 


91.4 


146.55/8 


17.74/5 


87.9 




2693/21 


1886/23 


30.0 


b 


3217/23 




Gerdt 3^= 


4596/21 


3051/23 


33.6 


5584/21 


4564/23 


18.3 


Arnborg-Lazard 


a 


a 




3042/15 


2476/11 


18.6 


Parametric Curve 


420.7/16 


17.72/10 


95.8 


522.5/16 


35.15/10 


93.3 


Katsura 4 


a 


a 




1059/13 


873/13 


21.1 


Arnold 1 


153.78/3 


a 




2276563/3 


264810/3 


88.4 


Arnold 2 


a 


a 




2919337/2 


1499222/2 


48.6 



a. Integer overflow. 

b. GMP problem. 

c. Garbage coflectioii turned off. 

d. Memory exhausted when garbage coUection is turned off. 

e. Total degree based ordering faster than prime-based ordering. 



With the "Katsura 4" example, we chanced upon a case in which the algorithm using 
total-degree ordering was faster than that with prime-based ordering. So, we have carried 

out a limited investigation of the effect of permuting the relative ordering of the indeter- 
minates of the two polynomial sets, for "Example 2" with three indeterminates, and for 
"Katsura 4" with five indeterminates. We found permutations in which the prime-based 
ordering was faster than the fastest total-degree ordering in both cases. 

The effects of permuting the relative order of the indeterminates for "Example 2" are 
presented in Tabic 2. In all cases using the 64 bit coefficients, the prime-based ordering 
version is faster by as much as 46%. When using GMP, the prime-based ordering is faster 
by as much as 30% in all but one case, when it is 16.7% slower. The fastest computation 
was the prime-based 64 bit case of 1.98ms when the indeterminate order was acb. The 
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ratio of maximum to minimum times is given for each implementation. For example, the 
ratio of the slowest prime-based 64 bit case to the fastest is 11. 



Table 2. Effect of the relative ordering of the indeterminates in "Exeimple 2" 



Order 


64-bit integer coefficients 


Multiple-precision integer coefficients 


Power Products 


Reduction 

(%) 


Power Products 


Reduction 

{%) 


total degree 


prime based 


total degree 


prime based 


abc 


15.85 


10.56 


33.6 


38.22 


33.57 


12.2 


acb 


2.93 


1.98 


32.4 


7.32 


6.4 


12.6 


bac 


24.06 


18.76 


22.0 


60.95 


58.73 


3.6 


bca 


41.48 


22.29 


46.3 


125.42 


88.22 


29.7 


cab 


3.71 


2.62 


29.4 


9.42 


8.36 


11.3 


cba 


15.08 


13.15 


12.8 


37.50 


45.01 


(16.7)" 


max/min 


14 


11 




17 


14 





a. Total degree based ordering faster than prime-based ordering. 



In the case of "Katsura 4" , there are five indeterminates, and 120 orders to consider. 
Only the multiple precision implementations work, as the 64-bit coefficient implementa- 
tions overflow. The fastest case is prime based, with a duration 17.6% less than fastest 
total degree based case. The ratio of the slowest to fastest is 1.9 for the total degree 
based implementation and 3 for the prime-based implementation. 

In the "Gerdt 3" example and some other examples, using the GMP library triggered 
a fault which caused the program to crash. This crash occurs when objects are collected 
by the garbage collector while running or when the program terminates. This is believed 
to be a problem in the interface between Eiffel and the GMP library rather than GPM 
itself. Turning ofi^ garbage collection avoided the fault, but proved costly, as it is faster 
to re-allocate memory internally than by repeatedly requesting additional memory from 
the operating system. Secondly, with garbage collection off, if memory is exhausted, the 
computation becomes disk-bound; these cases were abandoned. 

6. Conclusions 

Prime-based ordering, based on ordering power products by encoding the indetermi- 
nates as prime numbers and using the natural number ordc;r, is an admissible ordering. 
Prime-based ordering is not a lexicographical or total degree ordering. Implementations 
of this ordering reduce power product operations to integer operations. 

Several versions of Buchberger's improved algorithm have been developed and tested. 
Each result has been verified to satisfy the necessary and sufficient conditions to be 
a reduced Grobner Basis. Resulting bases also reduce their respective given polynomial 
sets, confirming that the Ideal is correctly preserved. They were also shown to be reduced 
by bases generated using the Grobner package in Maple. 
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Duration reductions measured using the improved Buchbcrger algorithm range from 
30% (40% speedup) to 96.8% (30x speedup). 

The number of polynomials can differ according to the ordering scheme. For example, 
in the "Gerdt 1" case, prime-based ordering generates a Grobner basis with 36 polyno- 
mials, whereas the total degree result has 56 polynomials. 

Finally, we have also explored the effect of permuting the indeterminates in some 
examples, and have found that the duration of computation varies significantly. In these 
examples, the fastest case was always using the prime-based ordering. We are continuing 
to investigate this matter together with issues associated with large coefficient size. 
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A. Polynomial Sets 

A.l. Example 1 

- 6ac3 + IZa'^c^ - I2a^c + Aa'^ 
6^ - 2ab - 2bc + + 2ac + 
+ 4:0 + 3 

ae" - ic^a^ + ica^ - + - ic^a + ica^ - + ac- 2a^ + 3c - 6a 
5abc^ + 3a&^ + a + 1 



A. 2. Example 2 

- &ac? + lia^c^ - \2a^c + Aa^ 
6^ - 2ah - 26c + + 2ac + (? 
+ Aa + Z 

ac" - Zc^a^ + 3ca3 - + - Zc^a + ica^ -a^ + ac- 2a^ 3c - 6a 
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A. 3. Example 3 



- 6ac^ + ISa^c^ - Ua^c + 4a^ 
+ 4o + 3 

ac^ - 3c^a^ + 3ca^ - a* + - 3c^a + 3ca^ -a^ + ac- 2a? + 3c - 6a 
hahc^ + 3a\? + a + 1 



A.J^. Cyclic 4 



X + y + z + t 
xy + yz + zt + tx 
xyz + yzt + ztx + txy 
xyzt — 1 



A. 5. Cyclic 5 



x+y+z+t+u 

xy + yz + zt + tu + ux 

xyz + yzt + ztu + tux + uxy 

xyzt + yztu + ztux + tuxy + uxyz 

xyztu — 1 



8 



A. 6. Gerdt 1 



yw — l/2zw + tw 

-2/7mu;2 + lO/lvw^ - + tu- 5tv + Otw 

2/7yw'^ " 2/7zw;2 + 6/7tw'^ - yt + zt - 
—2v^ + Auvw + bv'^w — 6uw'^ — Ivw^ + Ibw^ + A2yv 
-lAzv - 63yw + 21zw - A2tw + 147a; 

-9/7uw^ + A5/7vw^ - 135/7w;^ + 2zv'^ - 2tv^ - Azuw + lOtuw - 2zvw - 28tvw 
+4:zw'^ + 86tw'^ - A2yz + 142^ + A2yt - lAzt - 2lxu + 105.ot - 'il^xw 

%/7yw^ ~ %l7zvs- + ml7tvs- - 2xv^ - 4ytw + 6ztw - 24:t'^w + Axuw + 2xvw 
— 4xw^ + bQxy — ibxz + 84a:t 

2uvw — Qv^w — uw"^ + Vivvo^ — bw^ + lAyw — 2Stw 

v?w — 3uvw + buViP' + 14r/?j; — 28tw 

—2zuw — 2fMw + Ayvw + 6zvw — 2tvw — 16yw^ — lOzw^ + 22tw^ + A2xw 
28/3yuw + 8/'Szuu' — 20/3tuw — 88/3yvw — 8zvw + 68/3tvw + 52yw'^ + 40/3zw'^ 

-Mtw^ - 84:xw 
—4yzw + lOytw + 8ztw — 20t^w + 12xuw — 'iOxvw + IbxuP' 
—y'^w + l/2yzw + ytw — ztw + 2t^w — 3xuw + Gxvw — 3xw^ 
8xyw — Axzw + 8xtw 

A.l. Gerdt 2 

35y* - 30a;y2 - 210^^^; + 3x^ + 3^xz - 105z^ + UOyt - 21u 

5xy^ - UOy^z - 3x'^y + A5xyz - A20yz^ + 210yH - 25xt + 70zt + 126yu 

A. 8. Gerdt 3 

6xy^t ~ x^zt ~ Gxyzt + 3xz^t — 2z^t — 6xy^ + 6xyz — 2xz'^ 

-63xyH'^ + 9x'^zt^ + 63xyzt^ + I8y'^zt^ - 27xzH'^ - 18yzH'^ + l8zH^ + 78xyH 

-78xyzt - ISy^zt + 24xz'^t + \8yzH - 9zH - ISxj/^ + Ihxyz - hxz"^ 
18x^yH - 3x^zt - I8x'^yzt + Uxy^zt + 5x^zH - 12x^^2^ + 6xzH - 8zH - ISx^j/^ 

+18x^yz — 12xy^z — ix^z^ + 12xyz^ — (Sxz'^ 
-x^yt + 3xyH + lOyH - 15y'^zt + 3yz^t - 3xy^ - lOy^ + xyz + 15y^z - 5yz^ 
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A. 9. Arnborg-Lazard 

x^yz + xy^z + xyz"^ + xyz + xy + xz + yz 
x'^y'^z + x^yz + xy^z^ + xyz + x + yz + z 
x'^y'^z'^ + x^y^z + xy^z + xyz + xz + z + 1 

A. 10. Parametric Curve 

x^l — x^ — X — y 

x^ - z 
x^O-t 

A. 11. Katsura 4 

2x^ + 2?/2 + + + - u 
xy + 2yz + 2zt + 2tu-t 
2xz + 2yt + + 2zu- z 
2xt + 2zt + 2yu-y 
2x + 2y + 2z + 2t + u-l 

A. 12. Arnold 1 

Sx^y"^ + bxy^ + ix^z + x^yz 
x^ + 2y3z2 + Uy'^z'^ + ftyz^ 
8a;-'' + 12y3 + xz'^ + 3 
Tx^y^ + ISxy'^z^ + y^z'^ 

A. 13. Arnold 2 

2xy'^z^ + x^y'^z - x^y^z + 2xyz'^ + 7t/^ + 7 

2x'^y'^z + x'^yz'^ - xy^z'^ + 2x^yz - 12x + 12y 

2y^z + x'^y'^z - xy^z - xy^ + + 2y'^z 

Zxy^z'^ + x^y^z - xy^z + ^)^z'^ + 2>xyz^ + 42;^ -x^y 
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